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We study the transfer of spectral weight in the photoemission and optical spectra of strongly 
correlated electron systems. Within the LISA, that becomes exact in the limit of large lattice coor- 
dination, we consider and compare two models of correlated electrons, the Hubbard model and the 
periodic Anderson model. The results are discussed in regard of recent experiments. In the Hubbard 
model, we predict an anomalous enhancement optical spectral weight as a function of temperature 
in the correlated metallic state which is in qualitative agreement with optical measurements in V2O3. 
We argue that anomalies observed in the spectroscopy of the metal are connected to the proximity 
to a crossover region in the phase diagram of the model. In the insulating phase, we obtain an excel- 
lent agreement with the experimental data and present a detailed discussion on the role of magnetic 
frustration by studying the fc— resolved single particle spectra. The results for the periodic Anderson 
model are discussed in connection to recent experimental data of the Kondo insulators Ce^BiiPta 
and FeSi. The model can successfully explain the different energy scales that are associated to the 
thermal filling of the optical gap, which we also relate to corresponding changes in the density of 
states. The temperature dependence of the optical sum rule is obtained and its relevance for the 
interpretation of the experimental data discussed. Finally, we argue that the large scattering rate 
measured in Kondo insulators cannot be described by the periodic Anderson model. 

PACS numbers: 71.27.-ha, 71.28.-Hd, 78.20.Bh, 71.30. 4-h 
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I. INTRODUCTION. 

The interest in the distribution of spectral weight in 
the optical conductivity of correlated electron systems 
has been revived by the improvement ofithe quality of 
the experimental data in various systemsETcl. 

The traditional methods used in the strong correla- 
tion problem, exact diagonalization of small clustersQ, 
slave boson approachesO, and perturbative calculations, 
have not been very successful in describing the interesting 
transfer of optical weight which takes place as a function 
of temperature in the strong correlation regime. 

Recently, much progress has been achieved by mapping 
lattice models into impurity models embedded in an effec- 
tive medium. This technique, the Local Impurity Selfcon- 
sistent Approximation (LISA)eI, is a dynamical mean field 
theory that becom^ exact in the limit of large number of 
spatial dimensiongCl. For instance, the Hubbard and the 
Anderson lattice model can be mapped onto the Ander- 
son impurity model subject to different sejfponsistency 
conditions for the conduction electron bathEicl. These re- 
sulting selfconsistent impurity problems-can be analyzed 
by a variety of numerical techniquesEJ~t3. 

In this paper we apply this approach to the study of the 
optical conductivity in regard of the recent experiments 
in V2O3, Ce^Bi/^Ptz, and FeSi. We take the view that 
the low energy optical properties of V2O3 can be modeled 
by a one band Hubbard model, while Ce^Bi/^Etp, and 
FeSi are described by a periodic Anderson modelE3. The 
modeling of the experimental systems requires a large 
value of the Coulomb repulsion U. 

Our main goal in this work is to demonstrate that 
simplified models of strongly interacting systems treated 
with the LISA, can account for the main qualitative fea- 
tures that are observed experimentally in strongly cor- 
related electron compounds. The paper is organized as 
follows: in section ^ we summarize the mean field equa- 
tions for the model hamiltonians and the expressions for 
the calculation of the optical conductivity and the op- 
tical sum rule. In section HI we present and intuitive 



pedagogical discussion the physical content of the solu- 
tion of the model hamiltonians in the large dimensional 
limit. Section IV is dedicated to a thorough discussion 
of the optical conductivity results. We discuss the effects 
of magnetic frustration on the spectral functions of the 
Hubbard model. And the effects of temperature and dis- 
order on the optical spectra of the Anderson lattice. The 
theoretical calculations are carried out using exact diag- 
onalization (ED) and iterated perturbation theory (IPT) 
techniques, and compared with texperimental results on 
various systems. We stress that the use of IPT allow 
us to access physically interesting regimes which are out- 
side the scope of the Quantum Monte Carlo method, and 
that the exact diagonalization technique is used to con- 
firm that the results presented are genuine features of 
the Hubbard model in infinite dimensions and not arti- 
facts of the IPT. The conclusions are presented in the last 



section. Our ED approach to the solution of correlated 
models in large dimensions is based on the use continu- 
ous fractions. The Appendix describes an algorithm to 
convert the sum of two given continued fractions into a 
new continued fraction which we use to extend the ED 
method to the models we treat in this paper. Part of 
the theoreticaLxesults in section IV were announced in 
a recent letteiniJ . The optical conductivity of the An- 
derson model and the Hubbard model were considered 
previously by Jarrell et al. using the Daamtum Monte 
Carlo and Maximum Entropy methods .liJOS 



II. METHODOLOGY. 



A. Mean field equations. 

As model hamiltonians we consider the Hubbard and 
the periodic Anderson model (PAM): 

Hh = - Y + f^)cLcja + Y U{nti - hin^l - ^), 



<i,3> 



(1) 
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+ ^^l^- + + ^("'i'T - \){ndri - \) (2) 



where summation over repeated spin indices is assumed, 
/i is the chemical potential, and Uj is the hopping am- 
plitude between the conduction electron sites, which in 
the PAM results in the band e^- The S and d opera- 
tors create and destroy electrons on localized orbital with 
energy e^. is the hybridization amplitude between c 
and d— sites, which also appear in the literature as d and 
/—sites respectively. pp. 

The derivation has been given in detail elsewhereB'EI. 
So we only present the final expressions. The resulting 
local effective action reads, 

Siocai = - dT dr'^l{T)g^\T-T')Mr') 

Jo JQ 

+ U jyr{n^[r)-\){n,{T)~\) (3) 

where V'J, ^o- correspond to a particular site, and denote 
cj., Co- in the Hubbard model, and {cj.,(i^}, {cn,da} in 
the PAM case. Ua corresponds to rica and Uda respec- 
tively. Also note that Eq.|| defines the associated im- 
purity problem, with -i/jj., ip^r being the operators at the 
impurity site while the information on the hybridization 
with the environment is implicitly contained in Q^^ ■ Re- 
quiring that Giocaii'^) — '^kG{k, oj), we obtain as selfcon- 
sistency condition 
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g^\u;) =uj + ^l~t'^G{u;) (4) 
for the Hubbard model, and 

[go%,iLo)^uj + ti-t''[GUiio) (5) 
with Qq explicitly given by 

Go M=(^ 'y' ' (6) 

for the PAM. In both cases, G is the "cavity" Green 
function which has the information of the response of 
the lattice. 

We will consider the symmetric case with /i = and 
eJ5 = 0. Moreover, we assume a semi-circular bare 
density of states for the conduction electrons, p°{e) = 
•£Jie ~ ek)/N,ues - i2/7rD)^l - {e/Oy, with the 
half-bandwidth D = 2t. This density of states can be 
realized in a Bethe lattice and also ou_^ully connected 
fully frustrated version of the modeltJ'Ej. In this case 
the "cavity" Green function simply becomes G — G. In 
the following we set the half-bandwidth r£||-^ 1 ■ We use 
an exact diagonalization algorithm (ED)tZHl3 and an ex- 
tension of the second order iterative perturbatiori-theory 
(IPT) to solve the associated impurity problemtJ. We 
have checked that IPT and the ED method are in good 
agreement for all values of the model parameters. This 
results from the property of IPT ta capture the atomic 
limit exactly in the symmetric casellj. We use extensively 
the IPT on the real axis to scan through parameter space. 
A detailed comparison will be presented elsewhere. 



B. Optical conductivity. 

The optical conductivity of a given system is defined 

by 

1 r°° 

a{cj)^—lm / {[J{t),jm)e^-'dt (7) 

where V is the volume, is the current operator and ( ) 
indicates an average over a finite temperature ensamble 
or over the ground state at zero tempetaiijre. In general 
(t{uj) obeys a version of the f-sum rulecJ'Ej, 

aiiu)dLu = -lm{[P,J]) (8) 

where P is a polarization operator obeying ^ — J . 

In a model which includes all electrons and all bands 
the current operator is given by 

J = ^Y.p,5(r - n) (9) 

where pi is the momentum and the position of the i*''- 
electron, and e and m denote its charge and bare mass. 
P is given by 



P = eY,n5{r-n). (10) 

i 

Thus, -^{[PtJ]) = where n is the density of elec- 
trons, and the sum rule becomes 

ne^ 

I a(bj)dbJ = TT . (11) 

Jo rn 

This result is clearly temperature independent and does 
not depend either on the strength of the interactions. 

When dealing with strongly correlated electron sys- 
tems, in a frequency range where few of the bands are 
believed to be important, it is customary to work with an 
effective model with one or two bands, such as the Hub- 
bard or the periodic Anderson model. The current oper- 
ator is thus projected onto the low energy sector and is 
expressed in terms of creation and destruction operators 

of the relevant bands (i.e., J = i^J^ii^^l+s^i ~ ^l^i+s) 
for the Hubbard and Anderson model). It is in this case 

that the expectation value {[P,J']) is no longer ^ 
but, proportional to the expectation valutj of the kinetic 
energy (K) of the conduction electronsOc3. In general 
(K) depends on the temperature and strength of interac- 
tions, therefore, for these few bands models, the optical 
weight sum rule will depend on them as well. If the pro- 
jection onto a few band model is valid, this result also 
implicitly indicates that a portion of the optical spectral 
weight (the weight not exhausted by (K)) is transferred 
to much higher energies, that is, to the bands that were 
excluded by the projection to low energies. 

In this paper we do not address the question of the va- 
lidity of the low energy projection onto a few band model. 
Instead we focus on the consequences of this assumption 
on the redistribution of the optical weight within a mean 
field theory that is exact in the limit of large dimensions. 
Our main conclusion is that there is a considerable tem- 
perature dependence of the integrated spectral weight 
appearing in the sum rule. 

In infinite dimensions, (t{i-o) can be expressed in terms 
of the onp nartirlp spcctrum of the current carrying 
electronsEj'S 

<^i'^) = Ts— / de p {e) — 

A,{uj')A,{uj' + uj){nf{u;') -nf{uj' + uj) (12) 

with Ag{uj) — — 2Im[Gfe(a;)] being the spectral represen- 
tation of the Green function of the lattice conduction 
electrons, a the lattice constant, and v the volume of the 
unit cell. 

As we anticipated, the kinetic energy is related to the 
conductivity by the sum rule 

foc 2 2 2 

An important result, which will be demonstrated later 
on, is the notable dependence of the plasma frequency ujp 
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with temperature. This feature will be seen to emerge as 
consequence that correlation effects generate small en- 
ergy scales (e.g. the "Kondo temperature" of the associ- 
ated impurity). It is the competition between the small 
scales and the temperature that gives rise to an unusual 
temperature dependence to the integrated optical spec- 
tral weight. 

At r = 0, the optical conductivity of a metallic corre- 
lated electron system can be parametrized byEj 

UJ* ^ 

(t(w) = + areg {uj) (14) 

47r 

where the coefficient in front of the 5-function is the 
Drude weight and ujp is the renormalized plasma fre- 
quency. In the presence of disorder S{uj) is replaced by a 
lorentzian of width F. 

Evaluating equation (^2|) at T = one finds in mean 
field theory that, 

k 

where Z is the quasiparticle weight. For the Hubbard 
model in infinite dimensions the expression above further 
simplifies, and it only depends on the density of states 



III. PHYSICAL CONTENT OF THE MEAN 
FIELD THEORY. 

A. Hubbard Model. 

The solution of the mean field equations shows that at 
low temperatures the model has a metal insulator transi- 
tion (Mott-Hubbard transitipH)|-a|t an intermediate value 
of the interaction Uc ~ 3ii3il3. The metallic side is 
characterized by a density of states with a three peak 
structure: a central feature at the Fermi energy that 
narrows as one moves towards Uc from below, and two 
broader incoherent features that develop at ±-^, namely, 
the lower and upper Hubbard bands. They have a width 
« 2D and their spectral weight increases as the transition 
is approached. The insulator side, for U > Uc, presents 
only these last two high frequency features, which are 
separated by an excitation gap of size A w C/ — 2D. The 
different structures of the DOS (Fig.0) give rise to very 
different optical responses. 



Metal Insulator 




-U/2 +U/2 -U/2 +U/2 




U/2 U " U 

FIG. 1. Schematic DOS for the Hubbard model (1/2 filling) 
and their corresponding optical spectra for the metallic and 
insulator solutions. The width of the incoherent peaks in the 
DOS is ~ 2D and the one of central peak in the metal is 
^ ZD = e'p. 

Lets first consider the insulator, which is simpler. In 
this case, optical transitions are possible from the lower 
to the upper Hubbard band. We therefore expect the op- 
tical spectrum that results from the convolution (|l2|) to 
display a single broad feature that extends approximately 
from U - 2D to U + 2D (Fig.|l|). A neghgible temper- 
ature dependence of the spectra is expected, as long as 
T << A. On the other hand, in the metallic case, the 
low temperature optical spectrum displays various con- 
tributions: i) A narrow low frequency peak that is due 
to transitions within the quasiparticle resonance, in the 
T = limit this peaks becomes a 5— function and is the 
Drude part of the optical response, ii) At frequencies of 
order ^ a-n incoherent feature of width ~ 2D emerges due 
to transitions between the Hubbard bands and the cen- 
tral resonance. Hi) A last contribution at frequency ~ U 
appears due to transitions between the Hubbard bands. 
This is a broad feature of width ~ AD. Therefore, we 
expect an optical spectrum which is schematically drawn 
in Fig.|l|. It is important to realize that, unlike the insula- 
tor, a notable temperature dependence of the spectra is 
to be expected. There is a low energy scale Tcoh that 
corresponds to the temperature below which coherent 
quasiparticle excitations are sustained. It roughly cor- 
responds to the the width of the resonance at the Fermi 
energy ep = ZD. As T is then increased and becomes 
comparable to Tcoh, the quasiparticles are destroyed, and 
in consequence, the contributions to the optical spectra 
associated with them, (i) and (ii), rapidly decrease. 

It should be clear that in our previous discussion we 
have assumed that the system does not order magneti- 
cally, as paramagnetic solutions were considered. This 
situation can in fact be realized by introduction of dis- 
order (e.g. a random distribution of Uj) or next nearest 
neighbor hopping, and avoids-the artificial nesting prop- 
erty of the bipartite lattical^o. 
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B. Periodic Anderson Model. 

We now present a schematic discussion of the peri- 
odic Anderson model solution. In this case there are two 
different types of electrons, the c-electrons which form 
a band and the d-electrons with localized orbitals. In 
the non-interacting particle-hole symmetric case, the hy- 
bridization amplitude V opens a gap in the c— electron 
density of states Aind ^ V'^/D. On the other hand, 
the original (5— function peak of the localized d— electrons 
broadens by hybridizing with the conduction electrons 
and also opens a gap Aind- 

When the effect of the interaction term is considered, 
as the local repulsive U is increased, one finds that for low 
frequencies the non-interacting picture which was just de- 
scribed still holds, however, with the bare hybridization 
V being renormalized to a smaller value V*. Thus, we 
say we have a hybridization band insulator with the hy- 
bridization amplitude renormalized by interactions. This 
can also be interpreted by considering that the inter- 
acting d— electrons form a band of "Kondo-like" quasi- 
particles, that allows to define a coherence temperature 
T* similar to the Tcoh introduced before. This coher- 
ent band further opens a gap due to the periodicity of 
the lattice. This is the well known scenario that is born 
out from slp4e boson mean field theory and variational 
calculationgEa. On the other hand, the present dynami- 
cal mean field theory also captures the high energy part 
of the d— electron density of states that develops inco- 
herent satellite peaks at frequencies ±-j with spectral 
weight that is transferred from low frequencies. Conse- 
quently, the c— electron density of states is mainly made 
of a central broad band of half-width D — 2t and a gap 
at the center that gets narrower a,s V ^ V* . Also, it 
develops some small high frequency structures, that re- 
sult from the hybridization with the d-electrons. In Fig.|^ 
we schematically present the density of states for the c 
and d— electrons. As in the Hubbard model, we assume 
the absence of magnetic long range order (MLRO). For a 
study of the magnetic phase of the Anderson model see 
Ref. M 



Since the d— sites are localized orbitals, only the 
c— electrons contribute to the optical response of this sys- 
tem. At T = 0, following the previous interpretation 
in terms of a renormalized non-interacting hybridization 
band insulator and equation (]T^), we expect to find an 
optical conductivity spectra with a gap A^ir, which de- 
creases as the interaction is increased. We also expect 
that Aind << Adir, as the first corresponds to the indi- 
rect gap from the density of states Aind ^ V*'^/D, while 
the second is the direct gap A^ir ~ V* that is defined as 
the minimum energy for interband transitions at a given 
k (see Fig.|]). We do not expect any other important 
contributions to the optical response since, as we argued 
before, the incoherent high frequency structures of the 
c— electron density of states do not carry much spectral 
weight. In Fig.^we schematically present the optical re- 
sponse at T = 0. 

As the temperature is increased the gap in the opti- 
cal conductivity becomes gradually filled. At high tem- 
peratures a simple picture of electrons scattering ofi^ lo- 
cal moments emerge. The crossover between these two 
regimes, would naively occur at a temperature of the or- 
der of Aind- 

Thus, we note that in the Hubbard model and in the 
periodic Anderson model the destruction of a coherent 
quasiparticle state that sets the low energy scale of the 
system has rather opposite effects in the optical response. 
In the first case, the correlated metallic state is destroyed 
as T becomes of the order of the renormalized Fermi en- 
ergy, and the Drude part of the optical response is trans- 
ferred to higher energies as the insulating state sets in. 
In the second case, however, the destruction of the co- 
herent excitations is accompanied by the thermal closing 
of the gap in the density of states that turns the system 
metallic. As a consequence, the gap of the optical re- 
sponse is filled with spectral weight from higher energies 
to become a broad Drude-like feature. 



IV. RESULTS 



A. Hubbard model. 



p(co) 



-u/2 +U/2 

2D 



fV- 11.11 





a((0) 



" Adir k 

FIG. 2. Schematic DOS (1/2 filling) for c and d-electrons 
in the PAM (top). The corresponding schematic optical spec- 
tra at r = (bottom left) and the schematic band structure 
with the direct and indirect gaps (bottom right). 



We will discuss the results for the model in regard of 
different experimental data on the V2O3 system. Vana- 
dium oxide has three t2g orbitals per V atom which are 
filled with two electrons. Two electrons (one per V) are 
engaged in a strong cation-cation bond, leaving, the re- 
maining two in a twofold degenerate eg JaandEEI. LDA 
calculations give a bandwidth of - O.SeVEl. The Hub- 
bard model ignores the degeneracy of the band which is 
crucial in understanding the magnetic structureEJ, but 
captures the interplay of the electron-electron interac- 
tions and the kinetic energy. This delicate interplay of 
itinerancy and localization is responsible for many of the 
anomalous properties of this compound, and it is cor- 
rectly predicted by this simplified model. 
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Experimentally one can vary the parameters U and D, 
by introducing O and V vacancies or by applying presure 
or chemical substitution of the cation. We can use ex- 
perimental data to extract approximate parameters to be 
used as input to our model. In particular, from the exper- 
imental optical conductivity data in the insulating phase, 
a rather accurate determination can be made because, as 
it is apparent from the spectra, the lowj-frequency con- 
tribution is mainly due to a single peakEI. In regard of 
our schematic discussion of the previous section, the po- 
sition of the maximum should approximately correspond 
to the parameter U that corresponds to transitions from 
the lower to the upper Hubbard band. Also, according 
to the picture of the previous section, the total width of 
the peak should be ~ 4D which is twice the width of the 
Hubbard bands. Therefore, we can approximately esti- 
mate the parameter D as the distance from position of 
the peak maximum to the frequency where the feature 
decreased to half its height (see Fig.^. 




0.0 0.5 1.0 1.5 

M[eyi 

FIG. 3. The experimental (T(tj) of insulating V2-1/O3 with 
y = 0.013 at lOK (upper) and y = at 70K (lower). We 
indicate in the spectra the position of the maxima and their 
width from which the parameters U and D for the model 
calculations are extracted. 

The parameters from the metallic optical conductivity 
spectra are not so easily extracted. However, we can still 
obtain a rather precise determination by considering the 
difference spectra between the data at 170K and at 300 A' 
(see inset of Fig.^). 




0.0 0.5 1.0 1.5 

CO [eV] 

FIG. 4. The experimental o-{u)) of metallic V2O3 
at r = 170K (upper) and T = 300K (lower). 
The inset contains the difference of the two spectra 
Aa{u}) = ai70K{to) — crsooif (a;). Diamonds indicate the mea- 
sured dc conductivity adc- 

As we shall later discuss in detail, it turns out that 
the feature that appears in the difference spectra at a 
frequency « OAeV can be associated with the parameter 
U/2. This is also intuitively suggested by the schematic 
discussion of the previous section, as this feature corre- 
sponds to the enhancement of transitions from the lower 
Hubbard band to the central resonance at the Fermi 
level and from the resonance to the upper Hubbard band 
which are at a distance ~ U/2. The value for the param- 
eter D w OAeV in the metallic phase was determined 
by noting that: i) a priori there is no reason to expect 
that it should be much different than in the insulating 
phase (unlike the parameter U which could be modified 
by screening); ii) it is consistent with the recent LDA 
calculation that gives a half-width of « 0.5eV for the 
narrow bands at the Fermi levelEl; iii) despite the lack of 
very good experimental resolution the value is consistent 
with both the optical data that we reproduce in Fig.^and 
photoemission experimentsEj; iv) as will be shown later 
in the paper, this estimated value will allow to gather in 
a single semi-quantitative consistent picture the optical 
conductivity results with the V2O3 phase diagram and 
experimental results for the slope of the specific heat. 
The extracted parameters along with the values for the 
size of the optical gap (in the insulators) and the total 
optical spectral weight are summarized in table 



Phase 




Parameter 






D [eV] 


U [eV] A [eV] 


Lol/i-K [eV/ficm] 


Ins. (y=0) 


.33 ± .05 


1.3 ±.05 .64 ±.05 


170 ± 20 


Ins. (y=.013) 


.46 ± .05 


.98 ±.05 .08 ±.05 


800 ± 50 


Metal (170K) 


.4± .1 


.8±.l 


1700 ± 300 



TABLE I. Experimental parameters for the model. 



6 



In Fig]^ we display the phase diagram of the model in 
large dimensions generalized to include n.n. hopping ti 
and n.n.n. hopping t2- The condition t\+t^ = keeps 
the bare density of states p° invariant. For ^2/^1 — 
we recover the original hamiltonian, and t2/ti = 1 gives 
the paramagnetic solution. The cxtra^iopping provides 
a magnetically frustrating interactionts. This phase di- 
agram obtained for 12/ ti =j— »/l/3 has the same topol- 
ogy as the experimental oneE3t3. Frustration losvers the 
TNeei well below the second order Tmit pointE£|. Using 
the parameters of table |, Tmit ~ 24QK, which is only 
within less than a factor of 2 from the experimental re- 
sult. The dotted line indicates a crossover separating a 
good metal at low T and a semiconductor at higher T. 
Between these states pdc{T) has an anomalous rapid in- 
crease, as is shown in the results of Fig||. The reason for 
this feature can be traced to the thermal destruction of 
the coherent central quasiparticle peak in the DOS. We 
find the behavior of Pdc{T) to be in good agr|eement with 
the experimental results of Mc Whan et a].Ej. Another 
crossover is indicated by a shaded area, it separates a 
semiconducting region with a gap A comparable with T 
from a good insulator where T << A, consequently the 
crossover temperature increases linearly with U and the 
horizontal width of the crossover region becomes broader 
with increasing T . This crossover is characterized by a 
sudden increase in pdc as function of [/ at a fix T (in- 
set Figjsl). This crossover behavior wap experimentally 
observed in V2O3 by Kuwamoto et alE3. 
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FIG. 5. Approximate phase diag ram for the model with 
n.n. and n.n.n. hopping {12/ ti) = -^1/3. The 1st order para- 
magnetic metal-insulator transition ends at the critical point 
Tmit (square). The dotted line and the shaded region de- 
scribe two crossovers as discussed in the text. The full cir- 
cles indicate the position of the optical spectra. A: insulator 
{y = 0), B: insulator (y = .013), C: metal (y = 0, 170A"), 
D: metal [y = 0, 300A). Note that for comparison with 
experimental resuttp-increasing U/D is associated with de- 
creasing pressure.yy Inset; pdc{U) for T = 0.06D (full) and 
T = O.ISD (dotted). 




0.0 0.1 0.2 

T 

FIG. 6. Pdc{T) for U/D = 2.1,2.3,2.5 (bottom to top). 
The maxima of pdc{T) defines the dotted line. Obtained with 
the IPT method. 



1. Magnetically ordered solutions. 

In infinite dimensions the optical conductivity is a 
weighted convolution of two one particle spectral func- 
tions. The one particle spectral function is, therefore, 
the basic building block which gives rise to the various 
features of the optical conductivity. In this subsection we 
consider the nature of the spectral functions with mag- 
netic long range order (MLRO). The understanding of 
the qualitative differences and similarities between solu- 
tions with and without MLRO is relevant in regard of 
systems, like V2O3, that present both antiferromagnetic 
(AFI) and paramagnetic (PI) insulating phases. Results 
for the corresponding optical spectra will follow in the 
next section. 

In Figs and || we respectively show the single particle 
spectra of the PI and AFI insulating solutions for differ- 
ent values of the interaction U . The results are obtained 
from the ED method at T = for clusters of seven sites. 
The finite number of poles in the spectra correspond to 
the finite size of the clusters that can be practically con- 
sidered. A finite broadening of the poles was added for 
better visualization. 
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CO 

FIG. 7. DOS of the paramagnetic insulator solution ob- 
tained from exact diagonalization of 7 sites with U = 3, 5, 7 
(top to bottom). A small broadening has been added to the 
poles. 




CO 

FIG. 8. DOS of the antiferromagnetic insulator solution 
obtained from exact diagonalization of 7 sites with 17 = 3, 5, 7 
(top to bottom). A small broadening has been added to the 
poles. 



In the AFI case, we plot 
sublattice Green functions G, 



averaged value of the 



Ga 



1 



{GA(J + Gsa) — -^{GAa + Ga -ct) 



(17) 



which is the quantity to be compared to photoemission 
experiments. 

It is interesting to realize from these results, which cor- 
respond to rather large values of the interaction J7, that 
the spectra in both cases are roughly similar. They both 
present a lower and upper Hubbard bands at energies 
~ ±-j with a bandwidth « 2D and a corresponding gap 
A « C/ - 2D. 

In particular, the PI solution, merely presents a rigid 
shift of the incoherent Hubbard bands as the interaction 
U is varied, P^hi^h is reminiscent of Hubbard's solution 
to the modelEjO. On the other hand, in the AFI case, 
the shape of the density of states follows from the fact 
that the sublattice magnetization is basically saturated 



at these large values of the interaction. 

At [/ = 7, the largest value of the interactions con- 
sidered, we observe that the shape of the spectra of G 
becomes very similar to the corresponding one in the 
disordered case. This can be understood from the fact 
that the magnetic exchange scale J ^ ^ vanishes as 
U becomes large. As one decreases the strength of the 
interaction, we observe that the AFI spectra become in- 
creasingly different from the PI ones. In the former there 
is a transfer of spectral weight that occurs within the 
bands, from the higher to the lower frequencies. This is 
a consequence of the fact that, as the scale J becomes in- 
creasingly relevant, the spectra acquire a more coherent 
character. The "piling up" that occurs with the transfer 
of spectral weight as U is reduced, is the precursor of the 
weak coupling inverse square root singularity in the low 
frequency part of the density of states. It is interesting 
to note that recent photoemission experiments in V2O3 
report the presence of a small anomalous enhancement 
in the lower frequency edge of the spectrum in the AFI 
phase. This feature may be interpreted from the previ- 
ous results as evidence of the transfer of weight within 
the Hubbard bands. 

A complementary perspective on the results that we 
just discussed is obtained by looking at the fc-resolved 
spectra given by the imaginary part of the Green function 
G{k, uj) which reads. 



G{k,u) 



1 



(18) 



In the large limit the Green functions are labeled by 
the energy aO. Nevertheless, one can still think of this 
quantity as the analogous of the fc— resolved spectra if 
one notes that the e goes from —D to D as it traverses the 
band (the dispersion is linear in the non interacting case) . 
Thus, we can associate the "zone center vector" e = 
with the r point of the BZ, and the "nesting vectors" 
e = ±D with the commensurate M point. In Figs.^ and 
[lO| we show the e-resolved spectra for the different values 
of the interaction considered before. From the inspection 
of the spectra we observe that in the low J case for [/ = 7 
the single particle spectra remain basically unmodified as 
we scan the e "wave vectors" , which indicates the very 
incoherent character of the single particle excitations. On 
the other hand, as we lower U and the scale J becomes 
larger, we note the emergence of a peak in the U = 3 case 
for small values of e. Notice, also, the more dispersive 
character of the excitations. 
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FIG. 9. e-resolved single particle spectra of the antiferro- 
magnetic insulator solution obtained from exact diagonaliza- 
tion of 7 sites with U = 7 with e = 0.0, 0.25, 0.5, 0.75, 1.0 (top 
to bottom). A small broadening 77 = 0.1 was added to the 
poles and the figures were vertically shifted for better visual- 
ization. 
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FIG. 10. e-resolved single particle spectra of the antiferro- 
magnetic insulator solution obtained from exact diagonaliza- 
tion of 7 sites with U = 3 with e = 0.0, 0.25, 0.5, 0.75, 1.0 (top 
to bottom). A small broadening = 0.1 was added to the 
poles and the figures were vertically shifted for better visual- 
ization. 

An important parameter of the theory is the degree of 
magnetic frustration. Frustration is necessary_aot only 
to obtain the observed phase diagram of V2 03lij but to 
account for the general, shape of its angular integrated 
photoemission spectraES. We can summarize the results 
of this section by saying that the ED solutions indicate 
that as the degree of frustration is reduced and as U/t 
is reduced, the spectral function develops more disper- 
sion, and the excitations at low energy are more coher- 
ent (i.e. the imaginary part of the self energy is smaller). 
Many experiments place V2O3 in the regime of strong 
frustration, while the observation of dispersive features in 
the insulating phase of NiSi.^Seo.^SB may be explained 



by a lower degree of magnetic frustration in this com- 
pound. We shall now briefly consider the antiferromag- 
netic metallic state (AFM) that occurs upon the intro- 
duction of next nearest neighbor hopping in the nested 
lattice (i.e. partial frustration) at small (but finite) values 
of U (see phase diagram in Fig j|) . We obtain the density 
of states for is/ii = (AFI), and i2/ti = VT/S (AFM) 
with the interaction [/ = 1.5. The results obtained from 7 
sites exact diagonalization are shown in Fig.|ll|. It is very 
interesting to note that the peak structure of the density 
of states seems to be divided into low frequency features 
near cj = 0, and higher frequency structures at energies 
of the order of the band-width (which is also compara- 
ble to U for the chosen parameters). This is even more 
clear in the antiferromagnetic metallic state with partial 
frustration. 
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FIG. 11. DOS obtained from exact diagonalization of 7 
sites with U = 1.5 for (t2/ti) = (top) and (t2/fi) = \/l73 
(bottom). A small broadening 77 = 0.1 was added to the poles. 

We note that our results are qualitatively similar to the 
recent exact diagonalization results for the t ^ J model 
and also to quantum Monte Carlo results for the Hub- 
bard model on 2-dimensional finite size lattices with— a 
choice of parameters comparable to the one used herec^l. 
Furthemore, the similarity to the physics found in finite 
dimensional finite size lattices becomes more striking by 
comparing the e-resolved spectra of Fig.^ with the ones 
obtained by Preuss et al. in a recent QMC studycll. 
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FIG. 12. e-resolved single particle spectra of the antiferro- 
magnetic insulator solution obtained from exact diagonaliza- 
tion of 7 sites with U = 1.5 with e = 0.0,0.25,0.5,0.75,1.0 
(top to bottom). A small broadening rj = 0.2 was added 
to the poles and the figures were vertically shifted for better 
visualization. 

An important technical remark is that in order to ap- 
ply the exact diagonalization method of Ref. |l^ to the 
problem with intermediate frustration < t2/ti < 1, it 
is necessary to be able to average the continued fractions 
for the spin-up and spin-down Green functions into a sin- 
gle continued fraction. To perform this task we use the 
algorithm detailed in the Appendix. 

2. The insulating state. 

We now turn to the optical conductivity results. The 
experimental optical spectrum of the insulator was re- 
produced in Fig.Hl3. It is characterized by an excitation 
gap at low energies, followed by an incoherent feature 
that corresponds to charge excitations of mainly Vana- 
dium characterEI. These data are to be compared with 
the model results of Fig.^ The overall shape of the 
spectrum is found to be in very good agreement with the 
experimental results for the pure V2O3 sample. We dis- 
play the optical spectra results from both IPT and the 
ED method. The data show the very good agreement 
between this two methods. The peak structure in the 
ED data is due to the finite number of poles that result 
from the finite size of the clusters that can be considered 
in practice. 
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FIG. 13. The model (t(cij) for the insulating solution results 
at f/ = 4_D and T = from ED (thin) and IPT (bold). 

In Fig.|lJ we display the results for the size of the gaps 
A, which are in excellent agreement witij the experimen- 
tal results indicated by black squareai3. It is interest- 
ing to note that the results of Fig.|l^, shown for various 
degrees of magnetic frustration, indicate that in V2O3 
frustration plays an important role. The experimental 
system seems to be closer to the limit of strong frustra- 
tion, which is consistent with neutron scattering results 
that indicate different signs for thejinagnetic interactions 
between different neighboring siteso. 
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FIG. 14. The gap A versus U for the antiferromagnetic, 
partially frustrated and paramagnetic insulators (dotted, thin 
and bold). A is twice the energy of the lowest pole from the 
ED Green function. The data are for — » 00 from clusters 
of 3, 5 and 7 sites assuming l/ris scaling behavior. Black 
squares show the experimental gap for V2-y03 with y = 0.0 
and 0.013. 

Another interesting point is the fact that the gap ob- 
tained in the model optical spectra and the one obtained 
from the position of the poles in the single particle spec- 
tra coincide (Figs.p^and pT[). We therefore conclude that 
in this model the direct and indirect gap are very close 
(which justifies a posteriori that A is measured from the 
lowest pole of the local Green function). This result 
already predicted in Ref. |^, was experimentally can.- 
firmed by accurate recent photoemission study of V2O3L3. 
This follows from the fact that the imaginary part of the 
self-energy is very large wherever the electron density of 
states is non zero in the insulating solution (see Fig.p^. 
This is nothing but a direct consequence of the complete 
incoherent character of the upper and lower Hubbard 
bands. They describe a completely incoherent propaga- 
tion, and one should not think of them as usual metallic 
bands "shifted" by the interaction U . Notice that from 
the discussion in the previous section in the unfrustrated 
case, one expects a a larger difference between the direct 
and the indirect gap. 
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FIG. 15. lin[G{Lu)] and lm[T.{iu)] for [/ = 4 from IPT. Note 
that Im[E(Lj)] is large when Im[G(a;)] is non-zero indicating 
the incoherent character of the particle excitations. 

A final and important quantity that can be compared 

to the experiment is the integrated spectral weight 
which is related to (K) by the sum rule ([l^) . Setting the 
lattice constant a ~ 3A the average V — V distance, we 
find our results also in good agreement with the experi- 
ment (see Fig.nJ). 
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FIG. 16. Kinetic energy (K) versus f/ at T = for the 
antiferromagnetic and paramagnetic insulators (bold-dotted 
and thin), paramagnetic metal (bold), and partially frustrated 
model (thin-dotted). Black squares show the insulator exper- 
imental results. 

An interesting question, not yet fully settled, is the 
mechanism by which the insulating solution is destroyed. 
The destruction of the insulating state occurs at a point 
Uci which may be different from the critical point Uc ~ 
3D that is associated to the breakdorai-ef the metallic 
state as the interaction U is increasedE3ll3. This issue is 
physically relevant because one can envision a situation 
where the magnetic order stabilizes the insulating solu- 
tion over the metallic solution but due to a large degree 
of magnetic frustration, the insulating solution is very 
close to the fully frustrated paramagnetic insulator. The 
destruction of the paramagnetic insulating state was dis- 
cussed in Ref. |ll| using IPT. Here we address this issue 
using exact diagonalization. 



We first study the behavior of the gap in the one par- 
ticle excitation spectrum defined as the position of the 
lowest energy pole (with non negligible weight) in the 
Green function as a function of the number of sites in- 
cluded in the representation of the effective bath. Al- 
though the mean field theory is strictly formulated in 
the thermodynamic limit, in practice, the representation 
of the bath by a finite number of orbitals introduces fi- 
nite size effects. The data shown in Fig.|l^ were obtained 
from the extrapolation of results from finite size cluster 
Hamiltonians i?"" to the rig —> oo system. The value for 
A is defined as twice the energy of the lowest frequency 
pole appearing in the Green function. In Figjl^ we show 
the gap as a function of the interaction U in systems of 
Us = 3, 5 and 7 sites. Fig. 18 contains similar results as a 
function of l/us which shows the good scaling of A, es- 
pecially as the gap goes to zero as U is decreased. Thus, 
this approach indicates a continuous closure of the gap 
at a critical value of the interaction Ud ~ 2.15. 
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FIG. 17. The gap A versus the interaction U in the para- 
magnetic insulator. A is twice the energy of the lowest pole 
from the ED Green function. The data are from clusters of 
3, 5 and 7 sites (top to bottom). 
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FIG. 18. The gap A versus the inverse of the number of 
sites l/ris in the paramagnetic insulator for various values of 
U. A is twice the energy of the lowest pole from the ED Green 
function. 

We also investigate the behavior of the inverse mo- 
ments of the spectral function defined as: 



p{e)de 



(19) 



The behavior of these quantities give a more detailed 
picture of the transition. The local picture of the para- 
magnetic insulator is that of a spin embedded in an in- 
sulator. Hybridization with the bands of this insulator 
transfers spectral weight to high frequencies but the spin 
remains well defined at low energies (even though with a 
reduced spectral weight) as long as there is a finite gap 
in the insulator. As Uci is approached, and the gap de- 
creases we face the question whether the spin remains 
well defined even at the transition point. This depends 
on the behavior of the density of states of the bath pbath 
at low frequecies (we recall, pbath is essentiallyj-jo in a 
Bethe lattice, c/. Eq.|). Whittoff and Fradkiny have 
shown that if the density of states of the bath vanishes 
as a power-law pbath oc the spin remains well defined 
if /3 > 1 while the spin is Kondo quenched if /3 < 1 and 
the spin degree of freedom is absorbed by the conduction 
electrons. The case (3 = 1 is-marginal. 

In a previous publicationE3 we showed that within IPT 
the second inverse moment remains finite at the transi- 
tion, while it diverges in the Hubbard HI solution. 

Notice that m_2 can remain finite up to the transition 
even when the gap closes, but a divergent second inverse 
moment implies the continuous closure of the gap. In Fig. 
|T9| we plot the inverse of m_2 together with that of the 
first and third inverse moments. The results correspond 
to the extrapolation to the infinite size effective bath, 
performed similarly as was done previously for the gap. 
The inverse of the second inverse moment shows good 
scaling behavior with the system size and is found to 
go to zero for U « 2.12. At this value of the interaction 
the moment diverges, which signals the breakdown of the 
insulating state, with the gap closing continuously. As 
expected, the first inverse moment remains finite at the 
transition (it also shows good scaling behavior) and, on 
the other hand, the inverse of the third inverse moment 
becomes negative even before the transition. This is due 
to the fast divergence of the third moment which renders 
the finite size scaling inaccurate. It is important to stress 
that this way of looking at the transition is very different 
from the previous one, nevertheless, the estimates for Ud 
that are obtained after the infinite size bath extrapolation 
are consistently predicted to within less than 2%. The 
results are substantially different from the ones obtained 
from IPT. 
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FIG. 19. Inverse of the first three inverse moments 
{m-n)~^ of the density of states as a function of U. The 
three curves correspond from top to bottom to the inverse of 
the first, second and third inverse moment respectively. The 
results are the Ua ^ oo extrapolation from clusters of 3, 5 
and 7 sites assuming l/ris scaling behavior. The dotted con- 
tinuation of the last curve indicates results where the scaling 
is not reliable due to the strong divergent behavior of that 
inverse moment. The vertical dotted line indicates the value 
obtained for Uci from the inverse moment analysis. 



3. The metallic state. 

We now discuss the data in the metallic phase. In 
Fig.^ at the beginning of this section, we reproduced ex- 
perimental data for pwe samples that become insulat- 
ing below Tc ~ ISOi^EJ. The spectra were obtained in 
the metalhc phase at T = 170K and T = SOOif and 
are made up of broad absorption at higher frequencies 
and some phonon lines in the far infrared. They appear 
to be rather featureless, however, upon considering their 
difference (in which the phonons are approximately elim- 
inated) distinct features are observed. As T is lowered, 
there is an enhancement of the spectrum at intermediate 
frequencies of order O.beV; and more notably, a sharp low 
frequency feature emerges that extends from to O.lSe^. 
Moreover, these enhancements result in an anomalous 

change of the total spectral weight with T. We argue 
below, that these observations can be accounted by the 
Hubbard model treated in mean field theory. 

In Fig, 20 we show the calculated optical spectra ob- 
tained from IPT for two different values of T. The in- 
teraction is set to U = 2.1D that places the system in 
the correlated metallic state. It is clear that, at least, 
the qualitative aspect of the physics is already captured 
and setting D « OAeV we find these results consistent 
with the experimental data on V2O3 (Fig.^. As the tem- 
perature is lowered, we observe the enhancement of the 
incoherent structures at intermediate frequencies of the 
order ^■ to U and the rapid emergence of a feature at the 



lower end of the spectrum. This two emerging features 
can be interpreted from the qualitative picture that was 
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discussed in Sec. [11 which is relevant for low T. From the 



model calculations with the parameters of table I we find 
the enhancement of the spectral weight taking place at 
a scale Tcoh ~ 0.05D « 24QK which correlates well with 
the experimental data. Tcoh has the physical meaning of 
the temperature below which the Fermi liquid descrip- 
tion applies ,£3 as the quasiparticle resonance emerged in 
the density of states. 
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FIG. 20. The model cr(tj) for the metallic solution at 
U = 2. ID and T = O.OSD (upper) and 0.083D (lower). A 
small r — 0.3 and 0.5D was included to mimic a finite amount 
of disorder. 

In Fig.|l| we present the results for (K) as a func- 
tion of the temperature. An interesting prediction of the 
model is the anomalous increase of the integrated spec- 

tral weight ^ as T is decreased, a feature that is actually 
observed in the experimental data (note that the spec- 
tral weight is not recovered upto the highest frequencies 
where experimental data is available (w « 6eV). This 
effect is due to the rather strong T dependence of the 

kinetic energy (K) cx in the region near the crossover 
indicated by a dotted line in the phase diagram (Fig||). 
It results from the competition between two small en- 
ergy scales, namely, the temperature and the renormal- 
ized Fermi energy ep. 
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FIG. 21. Expectation value of minus the kinetic energy (K) 
as a function of the temperature for various U (IPT). This 
quantity is directly proportional to the optical conductivity 
sum rule. It predicts a notable increase in the optical spec- 
tral weight as the temperature is decreased in the correlated 
metallic regime. 

Fig.|2^ contains the comparison between the results for 
the same quantity (K) at f/ = 2 as obtained from the IPT 
and the finite temperature ED method. It demonstrates 
that the temperature dependence is indeed a true feature 
of the model which is being successfully captured by the 
approximate IPT solution. 

Although the qualitative aspect seems to be very ac- 

^2 

curately described by the model, we find -r^ « 10007=^ 
which is somewhat lower than the experimental result. 
This could be due to the contribution from tails of bands 
at higher energies that are not included in our model, or 
it may indicate that the bands near the Fermi level are 
degenerate. 
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FIG. 22. Comparison of the expectation value of the kinetic 
energy (K) as a function of the temperature for (7 = 2 as 
obtained from IPT and ED method. 

We now want to finally consider an important predic- 
tion of the model for the slope of the linear term in 
the specific heat 7 in the metallic phase. Experiments 
show that the slope 7 is in general unusually large. For 
0.08 Ti substitution 7 « 40 "^'^^ , while for a pressure 
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mj 



and with 



of 25Kbar in the pure compound 7 
V deficiency itua range of y 
7«47 

weight in the Drude peak in the optical conductivity and 

Theval- 



molK^ 

0.013 to 0.033 the value is 
^{^i -EJ In our model 7 is simply related to the 



— —3 

-7 n 



mJeV 



to the quasiparticle residue Z , 7 — ^ moiK'^ 
ues of [/ = 2. ID and D k, OAeV extracted from the op- 
tical data correspond to a quasiparticle residue Z w 0.3, 
and result in 7 w 25 ^^^ij which is close to the experi- 
mental findings. Thus, it turns out that the mean field 
theory of the Mott transition explains in a natural and 
qualitative manner, the experimentally observed optical 
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conductivity spectrum, the anomalously large values of 
the slope of the specific heat 7, and the dc-conductivity 
in the metallic phase, as consequence of the appearance 
of a single small energy scale, the renormalized Fermi 
energy e^. 



B. Periodic Anderson model. 

1. Gap formation. 

A second class of systems where the correlations in- 
duce an anomalous temperature dependence are the 
Kondo insulators. While the most qualitative physics 
of these systepifi is well understood, several features re- 
main puzzlingHj. The charge gap measured in optical 
conductivity is largePn than the spin gap As measured 
in neutron scattering. The transport gap At obtained 
from the activation energy in dc-transport measurements 
is smaller than Ac. The gap Ac begins to open at a 
characteristic temperature T* '--^ ^ and becomes fully 
developed at a much smaller temperature of the order of 
Also, the gap is temperature independent below T* . 
In CesBiiPts it is found that Ac « 450K, A^ « 250K, 
T* « lOOK and the optical gap is completely depleted 
only below « 25Kn. On the other hand, qualitatively 
similar results were reported for FeSi, with Ac « lOOOiir, 
T* « 200|fir and the gap becomes depleted between 20 
and lOOi^a 

The mean field theory accounts for all these observa- 
tions. The low energy behavior of the one particle Green 
functions of the model can be understood as that of a 
non interacting system where the interaction U reduces 
the hybridization from its bare value to a renormalized 
value V* which decreases as U increases. In consequence, 
the gap in the optical conductivity decreases by the effect 
of correlations. However, the lineshape remains approxi- 
mately invariant, and merely changed by a rescaling fac- 
tor respect to the response of the non-interacting model. 
This is demonstrated by the plot of the optical conduc- 
tivity a{uj) for different values of U shown in Fig.^. The 
optical gap Ac is given by the direct gap A^ir of the renor- 
malized band structure. These results were obtained by 
IPT at T = and we checked in various cases that the 
results are in excellent agreement with the ED method. 
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FIG. 23. The optical conductivity spectra of the periodic 
Anderson model for values of the interaction U = 0.5, 1, 2, 3 
(right to left), keeping the hybridization V = 0.25 fixed. The 
inset shows the gap from the optical spectra Ac « Adir 
and the indirect gap Ai„d from the local density of states 
for V — 0.6. The slopes of these curves indicate that 
oc Aind and V* oc Adir in the strong correlation re- 
gion. 

We now consider the behavior of (t^lj) with tempera- 
ture. Fig.^ shows the optical conductivity for different 
temperatures with the parameters U = 3 and V = 0.25 
fixed. The gap is essentially temperature independent. It 
begins to form at T* « 0.02 ~ and is fully depleted 
only at temperatures of the order of ^ . We thus observe 
that the mean field theory is able to capture the qualita- 
tive aspect of the experimental results that we summa- 
rized before. This basically consists in the individualiza- 
tion of 3 different energy scales: the largest corresponds 
to the gap of the optical spectra Ac ~ Adir, an interme- 
diate scale T* ~ where this gap starts to form and 
quasiparticle features start to appear in the DOS and, a 
third and smaller scale Amd ~ ^ , which corresponds to 
the temperature where the optical gap gets completely 
depleted. As demonstrated in Fig|2^ where we plot the 
results for the density of states, that smallest scale Aind 
also indicates the temperature below which the gap in 
the density of states opens, and, thus, can be associated 
to the gap measured in dc-transport experiments At. 
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FIG. 24. The optical conductivity for the Anderson model 
at r = 0.001 (bold), 0.005,0.01,0.02 (dotted), 0.03 (thin). 
The interaction U — 3 and V = 0.25. Inset: The same quan- 
tity at r = 0.001 (bold), 0.005, 0.01, 0.02 (dotted), 0.03 (thin) 
with lorentzian random site disorder of width F — 0.05. 
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FIG. 25. Low frequency part of the density of states for the 
d and c— electrons (top and bottom) obtained from IPT at T 
= 0.001, 0.005, 0.01, 0.02, 0.03 for [/ = 3 and V = 0.25 (top to 
bottom for electrons and bottom to top for c— electrons). 
Inset: The density of states in the full frequency range at 
T = 0.001. 

In order to make a meaningful comparison with the 
experimental data, we have added the effects of disorder 
by putting a lorentzian distributed random site energy on 
the conduction electron band with width F = 0.05. The 
results are displayed in the inset of Fig.^ and they show 
that the introduction of disorder makes the overall shape 
of the spectra in closer agreement with the experimental 
resultaJ'cl ( for a d iscussion of the scattering involved, see 
subsection I V B 2 ) . Also, we observe that increasing the 
disorder reduces the temperature T* . 

In the following, we briefly address the question of the 
integrated total spectral weight. It has been noticed that 
the experimental results in both, CesBiiPts and Fe 
seem to violate the sum rule for the spectral weightH 
However, this point has been recently questioned, at least 
for the FeSi compoundcS. In order to contribute to the 
proper interpretation of the experimental data, it is im- 
portant to compute the kinetic energy of our model at 
finite temperature, which is directly related to the sum 
rule of Eq.|l^. The results from IPT are presented in 
Fig.^ which shows the notable dependence of the ki- 
netic energy with temperature and interaction strength 
(we plot the negative of (K) which is the quantity that 
enters Eq.|l3|) . In Fig.|^, we plot similar results obtained 
with the ED algorithm which demonstrates that the be- 
havior captured by the IPT calculation is indeed a true 
feature of the model. 
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FIG. 26. Expectation value of minus the kinetic energy (K) 
as a function of the temperature for U = 0, 2, 3, 4 and V — 0.4 
(bottom to top) as obtained from (IPT). This quantity is di- 
rectly related to the optical conductivity sum rule. It predicts 
a notable decrease in the total optical spectral weight as the 
temperature is decreased in the range below the maxima. 



<K> 




-0.18 
0.00 



FIG. 27. Comparison of the expectation value of the kinetic 
energy (K) in the PAM as a function of the temperature for 
U = 2 and V = 0.4 as obtained from IPT and ED method 
(dashed line and full line). 

As we have previously discussed for the the Hubbard 
model case, the strong correlation effects that renders 
(K) a function of the temperature implies that if the 
PAM is the relevant model for the systems at low en- 
ergies, then the results predict the behavior of the in- 
tegrated optical weight within the low frequence range. 
Actually, experimental data which is inferred from the 
Kramers-Kronig transformation of reflectivity measure- 
ments, can only be reliably obtained within a limited low 
frequency range of the order of a few eV . The behavior 
of {K(T)) in Fig.^ is non-monotonic. As we increase 
the temperature from zero, we observe that initially the 
kinetic energy decreases. This is a consequence of the 
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electron delocalization since the system becomes a metal 
as the smaU gap in the density of states is filled. The ki- 
netic energy then goes through a minimum and starts to 
increase as the temperature is further increased. This is 
simply due to the thermal excitation of electrons within 
the single conduction band. Correlations now play an 
irrelevant role as the temperature is higher than the co- 
herence temperature T* . When we study the behavior for 
different values of the interaction U in Fig|2^, we observe 
that the position of the minima (maxima in this figure as 
— {K{T)) is plotted), becomes smaller as U is increased. 
This can be understood simply as a consequence of the 
rcnormalization of the hybridization amplitude V — > V* . 

In regard to the experimental situation in the Kondo 
insulators, which indicate the apparent violation of the 
optical sum rule, our results give a plausible qualitative 
explanation for the observed behavior. In fact, for exper- 
imental data obtained at temperatures smaller than the 
size of the gap Ac and restricted to a finite low frequency 
range (which is in fact the actual situation), the model 
predicts the apparent "disappearance" of spectral weight 
as the temperature is decreased. 

We should also point out that although this simple 
model accounts, rather successfully, for the various en- 
ergy scales, it fails to provide an accurate reproduction 
of the detailed experimental lineshape. A complete ex- 
planation of the experimental results may need the con- 
sideration of additional sources of scattering, as will be 
discussed in the next subsection. 

To finish our discussion on the gap formation in the 
periodic Anderson model we shall present the results for 
the size of the various gaps that are obtained from the 
correlation functions. 

The first study of the periodic Anderson model in large 
dimensions was caKried out by Jarrell et al. using quan- 
tum Monte CarloEa. Our spectral functions and density 
of states are in general agreement with the early work 
in the region where the QMC and exact diagonalization 
method can be compared. A noticeable qualitative dif- 
ference is that we find the spin gap to be slightly but 
strictly smaller than the indirect gap when U ^ 0. 

In Fig.|2| we show the local spin and charge correla- 
tion functions along with the optical conductivity which 
shows the qualitative agreement with the experimental 
data of Ref. ||. We also compare in the inset the di- 
rect optical gap Adir, the indirect gap Amd relevant for 
transport properties, and the spin gap Ag obtained from 
the spin-spin correlation function. We find that Adir 
is consistently larger than A^, and that As is somewhat 
smaller than Aind- As expected when U = 0, A^ = Aind, 
but as U increases becomes smaller than unity and 
approaches the value 1/2 at i7 « 2. 




FIG. 28. The local spin-spin (bold) and charge-charge 
(thin) susceptibility from 7 sites ED. The optical conductivity 
from IPT (dotted). The parameters are U = 1 and V = 0.2. 
The y-axis is in arbitrary units. Inset: The direct gap form 
IPT (upper dotted line), the indirect gap (lower dotted line) 
and the spin gap (solid line) from 8 sites ED. The hybridiza- 
tion is V = 0.2. 



2. The scattering rate. 

In the previous section we have stressed the qualita- 
tive success of the mean field theory of the periodic An- 
derson model in connection with the gap formation in 
Kondo insulators like FeSi and Ce^iBi^Pt^. In this sec- 
tion, however, we will show that this approach cannot 
account for the large scattering rate measured in these 
materials, if one does not include the effects of disorder 
in the model. This is very surprising and is an indica- 
tion of the limitations of the one band periodic Anderson 
model for modeling these systems. 

The optical conductivity of Kondo insulators is - ex- 
cept for the gap which forms at low temperatures - al- 
most constant over a large frequency range extending to 
several times the width of the gap. The correspond- 
ing value of cr(w) is quite similar for all of the materi- 
als (typically 3000 - 4000 (Slcrn)"^ ) and depends only 
weakly on temperature. The related scattering rate can 
be estimated (at T > T*) by simple Drude model ar- 
guments: at zero frequency, we have a — "^^"^ . Here, 
n = a~^, where a denotes the lattice constant, m can 

2 

be obtained from the kinetic energy ^ 



27r h 



D, where 

p « ■^^Tp- Assuming a « lO^^^^m, the equations yield 
h/r ^ (Tq ^10^ {ncm)~^ D . Thus the measured values for 
(To imply a scattering rate which is of the order of the 
bandwidth (i ~ D, assuming h = 1). This should be 
compared with the scattering rate found in normal met- 
als like copper, which is three orders of magnitude smaller 
(i ^ 10-313). 

Since all experiments on Kondo insulators (and also on 
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many Kondo metals) observe (above the gap) the same 
order of magnitude for cr(u;) one should expect that there 
is a common mechanism involved. It is reasonable to 
assume that the scattering of conduction electrons by 
the localized electrons in the periodic Anderson model 
provides an explanation. To address this question, we 
calculated the effective scattering rate. This quantity is 



D 



Hkl — 



determined by the effective c-electron self-energy S 



(20) 



where S^d is the self-energy of the localized d— electrons, 
which enters the formula for the optical conductivity (|l^ ) 
via A^^iuo) = — 2Im (e//), — . The imaginary 

part of Y!"cJ^\i^) measures the scattering involved. In 
Fig.^ we plotted this quantity for the particle hole sym- 
metric case, V = 0.251), U = SD, and T = 0.11?. Since 
in this section we are not interested in the gap forma- 
tion the temperature was chosen to be well above the 
point where the gap starts to open (T* « 0.0251?). For 
comparison, the calculation was done with both ED and 
IPT. 
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It is clear from the plot that the scattering rate is much 
smaller than the bandwidth D and gives rise to an opti- 
cal conductivity which is smaller than the experimentally 
observed value by two orders of magnitude. This result 
remains valid away from particle hole symmetry and for 
different choices of V and U . If one ignores the self consis- 
tency condition, one would expect, based on the theory 
of the Kondo impurity model that as the temperature 
is lowered the scattering rate should grow towards the 
unitary limit, D , this growth, which is expected at low 
frequencies and low temperatures is preempted, in the 
lattice by the formation of the hybridization gap. 

For ^ 1 and half filling, the periodic Anderson 
model can be transformed into a Kondo lattice model 
by a Schrieffer- Wolff Transformation: 



2 ^ 

<ij>,<T 



'^1(7 ''JO- 



(21) 



where J = Here, Si describes a spin at site i. For 

a cross-check, we also examined this hamiltonian using 
the exact diagonalization method. For J — ^ D, which 
corresponds to V — 0.25 D and U = 3D, we find at 
T = 0.1 D a scattering rate ^ ~ — ImEcc(w = 0) « 
7 X 10 ^'^-D. This is again much less than required to 
explain the experimental data. 

We believe that the failure observed here within the 
present approach is a general shortcoming of the periodic 
Anderson model. A more realistic description has to take 
several crystal filed split bands and this could increase the 
finite frequency optical absortion. 



V. CONCLUSIONS. 

In this paper we have illustrated how the LISA, that 
becomes exact in the limit of large dimensions, can be 
used to study the physics of systems where the local in- 
teractions are strong and play a major role. In particular 
we have demonstrated that the Hubbard and the periodic 
Anderson model treated within this dynamical mean field 
theory can account for the main features of the tempera- 
ture dependent transfer of spectral weight in the optical 
conductivity spectra. In the case of V2O3 we found that 
the theory is able to account semi-quantitatively for the 
conductivity results in both the metallic and insulating 
states. In the former case it can also account for the 
topology and energy scales of the experimental phase di- 
agram as well as for the unusually big values observed 
in the slope of the specific heat 7. In the latter case, 
the theory seems to provide further insights in the role of 
the magnetic frustration. In this regard, we have studied 
in detail the predictions of the model for photoemission 
spectra in the insulating phase with long range order and 
noted that the present mean field theory indeed captures 
many aspects of the behavior encountered in the numer- 
ical studies of the model in low dimensions. 

For the Kondo insulators, we have seen most of the 
qualitative features of the observed behavior of optical 
spectra with temperature being captured in detail by our 
model treated in mean field theory. In particular we iden- 
tified the different energy scales that describe the ther- 
mal filling of the optical gap and how they relate to the 
changes in the single particle spectra. However, we saw 
that the periodic Anderson model is not able to explain 
the high scattering rate measured in Kondo insulators. 

We presented results for the temperature dependence 
of the optical sum rule in the strongly correlated models. 
While in the Hubbard case, the results capture the qual- 
itative change of the total spectral weight with temper- 
ature observed in the V2O3 system, our quantitative re- 
sults on the PAM may be relevant for the resolution of the 
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"missing" spectral weight controversy in opticaLexperi- 
ments on the insulators CezBi^Pt^ and FeSoBo. From 
a broader perspective it has turned out to be very illumi- 
nating to realize how the emergence of a small "Kondo" 
energy scale which is a correlation effect results in an un- 
usual temperature dependence of the projected optical 
sum rule (c/. Eq.[T3|) in the Hubbard and the periodic 
Anderson model. We have shown that in the former case 
the optical weight increases when the temperature is re- 
duced and the system becomes more metallic, while in 
the latter it decreases, as a consequence of the system 
opening a gap and becoming insulating. 

We finally stress that our mean field approach can be 
easily adapted to incorporate more realistic band struc- 
ture density of states and more complicated unit cells. 
These extensions would allow for a more precise quanti- 
tative description of these interesting systems. 
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APPENDIX A: ADDITION OF TWO 
CONTINUED FRACTIONS. 

In this appendix we present an algorithm that allows 
to sum two continued fractions into a single one. This is 
necessary for the implementation of the ED method in 
models with magnetic frustration or disorder, where var- 
ious Green functions have to be averaged and the result 
has to be expressed as a new continued fraction. The 
details of the ED method can be found in Ref. 

In the ED method an effective cluster hamiltonian H"^' 
oius sites is diagonalized. At T = only the groundstatc 
\gs) and the groundstate energy Eq need to be obtained, 
and thiSjXan be efficiently done by the modified Lanczos 
method.cZl The local Green function G{uj) is then ob- 
tained as a continued fraction. Actually one needs to 
compute two continued fractions G^{ijj) and G^{u)), for 

< and for a; > respectively. 



G(w) = G>(w) + G<H 



{9s\c 



1 



1 



with 



c\gs) (Al) 



u! — a, 







1 >_ 



G<H 



(gslcfclgs) 



(A2) 



where c and are the operators associated with the 
local site of 7?"=. The parameters af^^ and bf^^ de- 
fine the continued fractions and are obtained from the 
following iterative procedure. 



(mm 
(/f-ii/r-i) 



(A3) 



where a =>, < and I/q) = c'lgs), = c\gs) and 

i/r+i>-^"i/n-<i/r>-fef i/r-i) 



(A4) 



and in the beginning we set 5q = 0. 

Thus, we observe that the basis defined by the vec- 
tors I/") gives a tri- diagonal representation of which 
contains the a^'s along the main diagonal and the ^/hi's 
along the diagonals next to the main one. In the follow- 
ing we drop the index a to simplify the notation. We will 
explicitly restore it in the final result. 

Let's now address the problem of our current inter- 
est. We assume that we have computed two Green 
functions G^(u;) where the index /i may label, for in- 
stance, a spin. Our task is to obtain a new contin- 
ued fraction representation of the average Green function 
G(w) = i(G|(w) -I- Gi{uj)). The more general case of a 
weighted average can be trivially generalized from the 
present case which we consider for simplicity. 

We first note that, from the Lanczos procedure, H"= 
has (different) tri-diagonal representations in the two 
sub-basis defined by \fi^) (we have dropped the >, < la- 
bel to simplify notation). The representation is basically 
a matrix that contains the parameters a' s along the main 
diagonal and the b's along the two sub-diagonals. 

The algorithm is as follows: one first diagonalizes the 
two tri-diagonal representations of by computing all 
the eigenvalues and eigenvectors. This is not numerically 
costly since the matrices are in tri-diagonal form and it 
may be done by standard methods. 

An important result that can be easily demonstrated 
is that the eigenvalues e'^^ of the tri-diagonal matrices are 
the poles of their corresponding Green functions G^{lo). 
Furthermore, one can also show that 



M 



u\2 



(A5) 



where v"^ are the first component of the M eigenvectors 
of the tri-diagonal matrices. 

Thus, from the definition of the Green function, one 
immediately recognizes that the vector {w^, w^, w*^} 
is nothing but cjjjgs) expressed in a basis where is 
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diagonal (which is a sub-basis of the given sector's Hilbert 
space) . 

The final step consists in writing the hamiltonian in 
the basis direct product of the two sub-basis, which, of 
course, will also be a diagonal representation of ; and 
then bring it to its tri-diagonal representation through 
the steps (A3), (A4) starting from the vector defined by 
(restoring the >, < label) 



\f>)^{c\+cl)\gs) 



ui e VI 



= { 



M „,M+1 „.M+2 



(A6) 



Thus, the newly determined af^s and bf^s that result 
from this last step are the parameters of the continued 
fraction representation of (uj) (the parameters for 
are obtained in a completely analogous manner). 
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